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ABSTRACT 

We use the Herschel-PCFLAS science demonstration data to investigate the star-formation 
properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminos- 
ity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly 
starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, 
we find that some individual objects are associated with high far-infrared luminosities. How- 
ever, the far-infrared properties of the radio-loud population are statistically indistinguishable 
from those of a comparison population of radio-quiet galaxies matched in redshift and ri- 
band absolute magnitude. There is thus no evidence that the host galaxies of these largely 
low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a popula- 
tion, have particularly unusual star-formation histories. Models in which the AGN activity in 
higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict 
a luminosity-dependent effect that is not seen in our data (which only span a limited range in 
radio luminosity) but which may well be detectable with the full Herschel- ATLAS dataset. 

Key words: galaxies: active - radio continuum: galaxies - infrared: galaxies 
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1 INTRODUCTION 

It is well known that the increase in the star formation density of 
the Universe with redshift (e.g. Madau et al. 1996) is paralleled by 
an increase in the luminosity density of quasars (e.g. Boyle & Ter- 
levich 1998), suggesting a link between galaxy assembly and ac- 
cretion onto massive black holes. It is much less obvious whether 
this link is a direct one: is star-formation activity physically as- 
sociated with AGN activity, and, if so, with what types of AGN 
activity is it associated? Although modelling suggests that the link 
may be a simple causal one, in the sense that mergers can trigger 
both AGN activity and star formation (e.g. Granato et al. 2004; di 
Matteo et al. 2005) this is a question that must be settled by direct 
observation of the AGN and star-formation properties of large sam- 
ples of galaxies. Far-infrared observations, probing dust heated by 
young stars, provide one of the best methods of measuring the star- 
formation rate, but in the past it has been hard to obtain statistically 
robust samples, although this has improved recently (e.g. Serjeant 
& Hatziminaoglou 2009; Serjeant et al. 2010; Bonfield et al. 2010). 

Radio-loud active galaxies are an important population in the 
study of the relationship between star formation and AGN activity. 
At the lowest radio luminosities (a population which is best stud- 
ied at z ~ 0) the radio galaxy population is dominated by objects 
for which there is no evidence at any waveband for any radiatively 
efficient AGN activity, setting aside non-thermal emission associ- 
ated with the nuclear jet (e.g. Hardcastle et al. 2009 and references 
therein). These objects have traditionally been called low-excitation 
radio galaxies (e.g. Hine & Longair 1979; Laing et al. 1994; Jack- 
son & Rawlings 1997), but the differences between them and their 
'high-excitation' counterparts are not only a matter of emission- 
line strength, but extend to optical (Chiaberge et al. 2002), X-ray 
(Hardcastle et al. 2006) and mid-IR (Ogle et al. 2006; Hardcastle et 
al. 2009). The low-excitation objects, where the AGN power output 
is primarily kinetic, clearly do not have the capability to regulate 
and eventually terminate star formation through coupling of their 
radiative output to cold gas. Instead the power of the AGN is put 
into the expansion of radio lobes, and therefore predominantly does 
work on the hot phase of the interstellar medium. At the highest ra- 
dio luminosities, though, radio-loud AGN activity is almost always 
associated with radiatively efficient accretion, and these objects can 
have strong effects on both the hot and cold gas in their environ- 
ments (cf. the 'radio mode'/'quasar mode' dichotomy of Croton et 
al. 2006; powerful radio-loud AGN are operating in both modes 
simultaneously). 

The relationship between radio-loud AGN and star formation 
might thus be expected to be a complex one, depending on radio lu- 
minosity, redshift and possibly radiative efficiency of the AGN, and 
indeed existing data present a picture of the relationship that seems 
to depend strongly on what AGN population (or individual object) 
is selected and which star-formation indicator is used. Hardcastle 
et al. (2007) argued that the nuclear differences between low- and 
high-excitation radio galaxies, suggesting a difference in accretion 
mode, might be explained by different sources of the accreting ma- 
terial, with the low-excitation sources powered by accretion (direct 
or otherwise) of the hot phase of the IGM, while high-excitation 
sources would be powered by cold gas supplied by mergers. This 
model is quantitatively viable in specific cases (Hardcastle et al. 
2007; Balmaverde et al. 2008) and is qualitatively supported by the 
observation that low-excitation sources of similar radio powers pre- 
fer richer environments (Hardcastle 2004; Tasse et al. 2008). How- 
ever, it also makes a prediction that high-excitation sources will be 
preferentially associated with gas-rich mergers and therefore star 



formation, which is borne out by observations both at low redshift 
(Baldi & Capetti 2008) and at z ~ 0.5 (Herbert et al. 2010). 

Far-infrared/sub-mm studies of star formation in samples of 
radio galaxies have so far concentrated on high-redshift objects, in 
which emission at long observed wavelengths (e.g. 850 /im, 1.2 
mm) corresponds to rest-frame wavelengths around the expected 
peak of thermal dust emission (e.g. Archibald et al. 2001; Reuland 
et al. 2004). Working at these high redshifts with the available flux- 
limited samples in the radio necessarily restricts these studies to 
the most powerful radio-loud AGN. However, the availability of 
observations at shorter far-IR wavelengths with the Herschel Space 
Observatory (Pilbratt et al. 2010) opens up the possibility of studies 
of very large populations of more nearby objects. In this paper, we 
present an analysis of all the radio-selected objects identified with 
galaxies detected by Herschel in the 14 deg 2 field acquired in the 
Science Demonstration Phase (SDP) of the Herschel Astrophysical 
Terahertz Large Area Survey (H-ATLAS: Eales et al. 2010). Al- 
though the full area of H-ATLAS (550 deg 2 ) will be required to 
allow us to investigate the complete range of dependence of host 
galaxy properties on radio luminosity and redshift, we show that 
the new data shed some light on the nature of the hosts of the nu- 
merically dominant low-luminosity radio galaxy population. 

Throughout the paper we use a concordance cosmology with 
H = 70 km s" 1 Mpc" 1 , Q m = 0.3 and fl A = 0.7. 



2 THE DATA 

The following images and catalogues were available to us: 

(i) PSF-convolved, background-subtracted images of the SDP 
field at the wavelengths of 250, 350 and 500 pirn provided by the 
Spectral and Photometric Imaging Receiver (SPIRE) instrument on 
Herschel (Griffin et al. 2010). The construction of these images 
for the H-ATLAS SDP data is described in detail by Pascale et al. 
(2010) and Rigby et al. (2010). We do not consider the PACS (Pho- 
todetector Array Camera and Spectrometer; see Ibar et al. 2010) 
data here as they are not deep enough for detections of more than a 
small fraction of our targets. The SDP field consists of two obser- 
vations over a 16 deg 2 region of the sky centred at RA = 9 h 5 m 30 s , 
Dec = 0°5'0". We restricted our analysis to the sub-region of the 
SDP field in which there was good data from both Herschel scans 
(hereafter the 'good' area of the SDP field); this covers 14.38 deg 2 . 

(ii) A catalogue of FIR sources detected in the SDP field, which 
includes any source detected at 5<r or better at any SPIRE wave- 
length (Rigby et al. 2010). 

(iii) Radio source catalogues and images from the NRAO VLA 
Sky Survey (NVSS) (Condon et al. 1998) and Faint Images of the 
Radio Sky at Twenty-one centimetres (FIRST: Becker, White & 
Helfand 1995) surveys. These cover the whole of the SDP field. 

(iv) Catalogues and images from the United Kingdom Infra-Red 
Telescope Deep Sky Survey - Large Area Survey (UKIDSS-LAS, 
hereafter LAS: Lawrence et al. 2007). The LAS covers only 92 per 
cent of the area of the SDP field. 

(v) Redshifts from the Galaxy And Mass Assembly survey 
(GAMA: Driver et al. 2009, 2010). GAMA is a deep spectroscopic 
survey with a limiting depths of tab < 19.4 mag, z < 18.2 and 
Kab < 17.6 in the SDP field; details of the target selection and 
priorities are given by Baldry et al. (2010) and Robotham et al. 
(2010). As described in more detail in Smith et al. (2010: here- 
after S10), the GAMA catalogue for this area contains 12626 new 
spectroscopic redshifts in addition to 1673 redshifts from previous 
surveys in the area. 
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(vi) A catalogue of galaxies in the SDP field with photometric 
redshifts based on the LAS and Sloan Digital Sky Survey Data Re- 
lease 7 (SDSS-DR7: Abazajian et al. 2009) photometric data, as 
described by S10. We filter this catalogue so as to require a K-band 
detection, to exclude any sources that are point-like in either the 
LAS or SDSS parent catalogues, and to impose the magnitude se- 
lection r < 22 mag. 

(vii) A list of identifications between these galaxies and detected 
sources in the H- ATLAS data, again as described by S 10. For these 
identifications S10 define a reliability 7? which is a measure of 
whether a single optical (r-band) source dominates the observed 
FIR emission; they suggest that only sources with R > 0.8 be used 
for this to be the case. Throughout the paper we consider all sources 
in the S10 catalogue, but distinguish in our analysis between 'reli- 
able' (R > 0.8) and unreliable identifications. 

(viii) A catalogue of optically selected quasars from the SDSS 
(Schneider et al. 2010) and 2dF-SDSS Luminous Red Galaxy and 
Quasar (2SLAQ: Croom et al. 2009) surveys in the SDP area, gen- 
erated by Bonfield et al. (2010). We use this purely as a comparison 
population: the reader is referred to Bonfield et al. (2010) for details 
of the selection of these objects. 



3 THE RADIO-SELECTED SAMPLE 

We began by constructing a sample of radio-detected objects in the 
H-ATLAS SDP field. This was done in the following way: 

(i) We selected all catalogued NVSS sources in the 'good' area 
of the SDP field (796 in total). As the flux cutoff for the NVSS cat- 
alogue is 5cr (2.5 mjy), these are all clearly detected radio sources. 
The good short-baseline coverage of the NVSS data ensure that the 
NVSS flux densities are good estimates of the true total flux density 
of our targets. 

(ii) We then cross-matched to the LAS K-band images by over- 
laying radio contours on LAS images, accepting only sources 
which had an association between the FIRST or, in a very few 
(5) cases, NVSS radio images and a K-band object with the ap- 
pearance of a galaxy or a quasar. FIRST is used in preference to 
NVSS for identifications because of its much higher angular reso- 
lution, allowing less ambiguous identifications where compact ra- 
dio components are present. This process excludes some weak or 
diffuse NVSS sources where FIRST detections were not available 
and where the NVSS position is inadequate to allow an identifi- 
cation with a LAS source. In addition, where NVSS sources were 
found to be blends of two or more FIRST sources, we corrected the 
NVSS flux density by scaling it by the ratio of the FIRST flux of 
the nearest source to the total FIRST flux. As mentioned above, the 
LAS covers only 92% of the area of the SDP field, so the choice 
to use this as our reference catalogue slightly reduces our coverage 
but does not affect the sample completeness in any way. We have 
391 objects with LAS identifications. 

(iii) Finally, we cross-correlated with the catalogue of galaxies 
with photometric redshifts (S10) on the basis of the LAS associa- 
tion. This gave us a total of 187 radio-loud sources in the SDP field. 
Where a spectroscopic redshift was available from the GAMA cata- 
logue (including pre-existing SDSS redshifts) we used that in pref- 
erence to a photometric redshift in subsequent analysis. 94 of our 
sources (50 per cent) had spectroscopic redshifts determined in this 
way. The median error on the photometric redshifts for all objects 
is 0.03 (see S10 for a discussion of the errors). 



This process gives us a catalogue which is flux-limited in the 
radio (by virtue of the original selection from the NVSS) and also 
magnitude-limited in the optical (we require a K-band identifica- 
tion and also require r < 22). In practice, this means that we have 
few sources with z > 0.7 and none with z > 0.85. The catalogue 
is also likely to be strongly biased against radio-loud quasars since 
we have excluded sources that appear point-like from our galaxy 
catalogue. This has the desirable effect that the measured luminosi- 
ties will tend not to be strongly affected by beaming and that any 
contamination of the fluxes measured at Herschel wavelengths by 
non-thermal emission might be expected to be limited (cf. the re- 
sults of Hes et al. 1995). Analysis of the full catalogue without the 
restriction of identification with (relatively) bright optical galaxies 
will be presented elsewhere (Virdee et al. in prep.): here we concen- 
trate on the implications of the Herschel properties of these sources 
for the nature of radio-loud AGN activity. 

Fig. shows the radio luminosity-redshift plot for our radio- 
loud sample (here and throughout the paper we adopt a = 0.8, 
where S oc v~ a , for the JC-correction in the radio luminosity cal- 
culations; a = 0.8 is a typical observed value for low-frequency 
selected objects^], and we expect that the selection against point- 
like optical objects will tend to select against flat-spectrum radio 
sources). It will be seen that we probe a wide range of radio lumi- 
nosities. At the low-luminosity, low-redshift end, we expect from 
existing analysis of the local 1.4-GHz luminosity function (e.g. 
Mauch & Sadler 2007) that the population will be dominated by 
luminous star-forming galaxies rather than radio-loud AGN, al- 
though a few AGN may still be present. The starburst luminosity 
function cuts off steeply above a few x 10 23 W Hz -1 at 1.4 GHz, 
so we expect that almost all objects above this luminosity will be 
radio-loud AGN, since we have no reason in the parent catalogue to 
be biased towards starburst galaxies. The Fanaroff-Riley break, i.e. 
the luminosity at which the population of radio galaxies switches 
from being dominated by objects of Fanaroff & Riley (1974)'s mor- 
phological class I to mostly containing objects of Fanaroff-Riley 
class II (hereafter these classes are abbreviated FRI and FRII) is at 
1.2 x 10 25 W Hz" 1 at 1.4 GHz, so our sample will be numerically 
dominated by FRIs, i.e. objects that are traditionally thought of as 
low-luminosity radio galaxies, but still has a significant number of 
objects at FRII luminosities. We have made no attempt to classify 
the objects morphologically in the radio using the NVSS or FIRST 
data, as these datasets tend to lack the surface-brightness sensitivity 
needed for reliable classification; morphological investigations will 
be discussed in a future paper. 



4 THE FAR-IR PROPERTIES OF THE SAMPLE 

In this section we describe the properties of the sample in the far- 
IR. Throughout this section, FIR fluxes in the SPIRE bands are 
measured directly from the background-subtracted, PSF-convolved 
H-ATLAS SDP images described in Section|2] taking the flux den- 
sity to be the value in the image at the pixel corresponding most 
closely to the LAS position of our targets, with errors estimated 
from the corresponding noise map. As discussed by Pascale et al. 
(2010), PSF-convolved maps provide the maximum-likelihood es- 
timator for the flux density of an single isolated point source at a 



1 For example, the mean 178-750 MHz spectral index for the 3CRR sample 
is 0.79; see http://3crr.extragalactic.info/ . At these redshifts the calculation 
is insensitive to the exact value adopted. 
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Figure 1. Radio luminosity of the radio-loud sample as a function of red- 
shift. Sources nominally detected by Herschel (i.e. sources identified with 
Herschel objects in the S10 catalogue) are marked as filled green circles. 
If the association with the LAS galaxy is deemed 'reliable' in the S10 cat- 
alogue, the object is also marked with a red open circle. Any source not 
listed as detected in the catalogue, but detected down to the 2cr level in the 
250-fim images, is shown as a light blue star. Non-detections are marked 
with blue open squares. The solid line corresponds to a nominal 5<r flux 
cutoff in the NVSS of 2.5 mJy. Sources lying significantly below this line 
do so as a result of the deblending process described in the text. Dashed 
horizontal lines show radio luminosities corresponding to (from bottom to 
top) the dividing line adopted in the text between starbursts and AGN; the 
expected luminosity for a maximal (4000 Mq year - ) starburst; and the 
luminosity corresponding to the Fanaroff-Riley break. 

given position in the presence of thermal noise; this remains a rea- 
sonable approximation if there are small correlations between the 
positions of multiple sources, as we expect in real data due to phys- 
ical clustering of objects. The flux densities we measure slightly 
underestimate the total flux density if the source is resolved, as we 
have verified by comparing our flux densities with those in the cat- 
alogue of Rigby et al. (2010), but this is only likely to be a problem 
at the lowest redshifts. We make an approximate correction for the 
mean background due to confusion by subtracting the mean flux 
level of the whole PSF-convolved map from each flux density mea- 
surement. 

Throughout this section we distinguish between sources cat- 
alogued by Rigby et al. (2010) and lower-significance detections, 
down to the 2a level, determined by the ratio of the flux density 
and the error measured by us from the maps at 250 /im: this dis- 
tinction is only for the sake of illustration, in that it allows us to 
show the properties of sources weaker than those in the source cat- 
alogue, and the presence or absence of a 2a detection forms no part 
of the quantitative analysis in the paper. 

To convert between measured flux or luminosity densities at 
Herschel SPIRE wavelengths and total flux or luminosity in the 
far-IR band, as used in the literature, we have to adopt a model for 
the far-IR spectral energy distribution (SED), since in general we 
do not have enough data to fit models to each object. We investi- 
gated several possible SED templates, including an optically thick 



grey-body model whose parameters were determined from a fit to 
Arp 220 (T = 61. 7K, = 1.54, as used by Stevens et al. 2010), 
an optically thin grey body with parameters fitted to M82 (T = 44 
K, p — 1.55) and the model fitted to normal galaxies in the SDP 
field (the S10 sample) by Dye et al. (2010) (T = 26 K, P = 1.5; 
cf. the very similar results, on a somewhat different sample of H- 
ATLAS sources, obtained by Amblard et al. 2010). These different 
choices of model parameters can make a significant difference (up 
to an order of magnitude) to the inferred total fluxes or luminosi- 
ties, and also, via the if-correction, to their inferred dependence 
on redshift. Because of this, we chose to use the model adopted by 
Dye et al., which has the merit of being derived from a dataset that 
has considerable overlap with ours; for our current purposes the 
absolute normalization of the FIR luminosity is less important than 
the relative normalization, since our main conclusions will come 
from a comparison of the radio-loud objects with other samples. 
Detailed model fitting to the SEDs of radio-loud objects, and thus 
the possibility of considering such factors as evolution in tempera- 
ture with luminosity or redshift will have to await a larger sample 
with a higher detection rate and will probably require additional 
constraints, such as the use of the PACS data; we plan to address 
this in a future paper. For the present work, we integrate the Dye et 
al. model between 8 and 1000 /im to obtain the total FIR fluxes and 
luminosities, for consistency with the approach used by other Her- 
schel papers, bearing in mind that this integration almost certainly 
underestimates the total IR flux since it includes no component that 
radiates in the mid-IR. 

4.1 Detections 

As can be seen from Fig. Q] only 31 of the galaxies identified 
with our 187 sample objects have identifications with Herschel 
sources from the catalogue of S10; all but 6 of these are classed 
as 'reliable' identifications, as discussed in Section [2] The detec- 
tion fraction is high at low radio luminosities (20/26 sources be- 
low L1.4 = 5 x 10 23 W Hz" 1 are detected) but then drops off 
rapidly with increasing radio luminosity. Given the known be- 
haviour of the starburst and AGN luminosity functions (Section 
2) this is best explained in terms of a dominant starburst popu- 
lation at low radio luminosities in our radio-selected sample. To 
illustrate this, we have computed the parameter q, as originally 
defined by Helou et al. (1985), for each of the detected objects. 
q — log 10 (S'iR/(3.75 x 10 12 x S1.4), where Sir is the integrated 
far-IR flux in W m -2 and Si. 4 is in W Hz -1 m -2 . We estimate 
q from the 250-/im flux density using the grey-body model dis- 
cussed above, taking proper account of A'-correction. q is plotted 
against radio luminosity in Fig. [2] For star-forming galaxies, q by 
this definition is expected to have a typical value around 2.4 inde- 
pendent of radio luminosity, as found for example using Herschel 
data by Ivison et al. (2010) and Jarvis et al. (2010), as a result of 
the well-known radio-FIR correlation (e.g. Condon et al. 1991). 
We observe that most of the low-luminosity sources have q values 
in the range 2-3. However, most of the detected high-luminosity 
sources have values of q in the range 0-1.5, implying radio flux 
densities 1-2 orders of magnitude above the values expected from 
their FIR fluxes on the radio-FIR correlation, and the upper limits 
for non-detections are at a comparable level. We conclude that we 
are indeed seeing predominantly starbursts at L1.4 < 5 x 10 23 W 
Hz -1 but that above this we are predominantly detecting genuine 
radio-loud AGN. These results are unaltered if we use values of q 
determined from fits to the 70-500 /im SEDs of the detected Her- 
schel sources (Jarvis et al. 2010) and the values of q obtained are 
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Figure 2. The parameter q (as defined in the text) as a function of radio 
luminosity for the sample, q is estimated from the flux in the 25Q-fj,m maps 
in all cases as described in the text. Colours and symbols are as for Fig. 
[T] The dashed horizontal line gives the median q (2.25) for the 20 reliably 
detected sources with L1.4 < 5 X 10 23 W Hz -1 . A 2cr upper limit on q is 
plotted for undetected sources. 

also consistent with that more detailed analysis. In what follows 
we use L1.4 < 5 x 10 23 W Hz -1 as an approximate luminosity 
cutoff to exclude the bulk of objects with q values consistent with 
being standard starbursts. While we recognise that more luminous 
starbursts than this can and do exist, the small number of objects 
observed to have L1.4 > 5 x 10 23 W Hz -1 and q > 2 suggests 
that few are present in our sample. 

4.2 Stacking analysis 

Since the vast majority of our radio sources are undetected at the 
5<t limit of the H-ATLAS source catalogue (Rigby et al. 2010), we 
need to use statistical methods to calculate the properties of the 
source population. We elected to stack the sample in bins corre- 
sponding to redshift and radio luminosity to investigate the way in 
which the far-IR properties varied with those two parameters. 

To establish quantitatively whether sources in the bins were 
significantly detected, we measured flux densities from 100,000 
randomly chosen positions in the field; a Kolmogorov-Smirnov (K- 
S) test could then be used to see whether the sources from our sam- 
ple were consistent with being drawn from a population defined 
by the random positions. Using a K-S test rather than simply con- 
sidering the calculated uncertainties on the measured fluxes allows 
us to account for the non-Gaussian nature of the noise as a result 
of confusion. We chose the boundaries of our luminosity bins in 
such a way that the lowest-luminosity bin contained all sources 
with L < 5 x 10 23 W Hz -1 and we also placed a bin boundary at 
the nominal Fanaroff-Riley transition luminosity, while otherwise 
trying to keep roughly similar numbers of sources per bin. Results 
are given in Tables Q] and [2] which also show the mean flux den- 
sity at each wavelength for every bin. We see that the sources in 



the lowest luminosity bin and the lowest two redshift bins are very 
strongly distinguished from the random background population, as 
expected from Fig. Q] sources at intermediate radio luminosity or 
redshift are not distinguished from the background, but then sev- 
eral higher-luminosity or higher- z bins are distinguished from the 
background at significance levels ranging from 2 to > 3er, at least 
at 250 /im. As expected (since the beam is larger and the potential 
for confusion greater), the significance drops off with increasing 
wavelength. Importantly, these results still hold, though obviously 
with somewhat reduced significance, if we exclude the formally 
detected sources (those with identifications in the catalogue) from 
the analysis, as shown in the bottom halves of Tables Q] and [2] Set- 
ting aside the lowest-luminosity and lowest-redshift bins (z < 0.2, 
L1.4 < 5 x 10 23 W Hz -1 ), where the exclusion of detections re- 
moves almost all the sources, all the bins that are significantly de- 
tected with the inclusion of the formally detected sources are still 
detected at 95 per cent confidence or better on the K-S test at 250 
fim if those sources are excluded, and many bins are more signif- 
icantly detected than that; so at least at high radio luminosities or 
redshifts we are not simply seeing the effect of a very few far-IR- 
bright objects. We emphasise, though, that it is the K-S statistics 
that include the detected objects (left-hand side of these tables) that 
determine whether a given bin is actually detected. 

We were also able to compare with the properties of galaxies 
that were not identified with radio sources. To do this we selected 
all objects identified as galaxies in the photometric redshift cat- 
alogue of S10 that had a LAS detection and a determined spec- 
troscopic (GAMA) or photometric redshift, imposing the cutoff 
z < 0.85, but were not identified with radio sources and that lay on 
the 'good' SDP field area so that photometry in the Herschel bands 
was possible. This gave 59,817 galaxies in total. The population of 
radio-selected sources is significantly (> 99.5 per cent confidence) 
different from this general galaxy population as a whole on a K-S 
test at all SPIRE wavebands. However, if we break both populations 
down by redshift, using the same binning scheme as previously and 
imposing the additional requirement that the comparison galaxies 
lie in the same K-band absolute magnitude range as the radio-loud 
hosts (which restricts us to 38,700 galaxies), we see that this ef- 
fect is completely dominated by the sources with z < 0.2, i.e. 
the sources that we identify as star-forming galaxies (Table[3}- The 
higher-redshift objects, which we expect to be radio-loud AGN, are 
indistinguishable, statistically, in their FIR flux distribution from 
galaxies selected in a similar way but without radio counterparts, 
and their mean flux densities at each wavelength are also very sim- 
ilar. These similarities are understandable, since we know that, at 
least at the low radio luminosities we are considering, radio galaxy 
hosts very often have the appearance of passively evolving ellip- 
ticals with no evidence either for strong star formation or for an 
obscured, radiatively efficient AGN, and the K-band selection and 
magnitude cutoff of our comparison population will have the effect 
of including many similar objects. We return to this point below. 

4.3 Far-IR luminosity and star formation 

We estimated far-IR luminosities in our bins using a modified ver- 
sion of the method of Serjeant & Hatziminaoglou (2009). The value 
of this approach is that it aims to account for the intrinsic scatter 
in the flux density of objects in bins as well as the scatter due to 
(thermal and confusion) noise when the luminosities are averaged; 
the reader is referred to Serjeant & Hatziminaoglou (2009) for a 
detailed description of the method. For each redshift or luminosity 
bin, we characterised the intrinsic scatter in the measured flux den- 
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Table 1. Mean bin flux densities and K-S probabilities that the Herschel fluxes of objects in redshift bins are drawn from the background distribution, as 
a function of wavelength. Low probabilities (below 1 per cent) imply significant differences between the bin being considered and the distribution of flux 
densities measured from randomly selected positions in the sky, as described in the text. 



Catalogued 
sources? 


z ram 


»e Objects 
in bin 


Mean bin flux density (mjy) 
250 fim 350 /im 500 fim 


K-S probability (%) 
250 fim 350 /im 500 fim 


Included 


0.00 - 


-0.10 


15 


396.1 ± 1.7 


160.6 ± 1.9 


58.8 ±2.3 


< 10" 3 


< 10~ 3 


< 10" 3 




0.10- 


-0.20 


14 


125.1 ± 1.7 


52.2 ± 1.9 


18.0 ±2.4 


< 10" 3 


0.004 


0.7 




0.20- 


-0.35 


39 


13.3 ± 1.0 


7.2 ± 1.2 


2.7 ±1.4 


2.4 


2.6 


20.7 




0.35- 


-0.50 


40 


6.9 ±1.0 


5.6 ± 1.1 


2.2 ±1.4 


0.07 


2.3 


40.1 




0.50- 


-0.65 


43 


5.9 ±1.0 


4.7 ± 1.1 


3.6 ±1.4 


1.0 


1.3 


1.9 




0.65- 


-0.85 


36 


10.0 ±1.1 


5.5 ± 1.2 


4.3 ±1.5 


0.007 


3.9 


5.2 


Excluded 


0.00 - 


-0.10 


4 


133.7 ±3.3 


69.2 ±3.6 


26.3 ±4.5 


30.6 


49.1 


92.9 




0.10- 


-0.20 


6 


5.3 ±2.6 


6.0 ±3.0 


1.9 ±3.7 


14.9 


25.7 


33.1 




0.20- 


- 0.35 


36 


3.1 ± 1.1 


3.6 ± 1.2 


1.9 ±1.5 


14.9 


18.8 


30.5 




0.35- 


-0.50 


37 


4.1 ± 1.1 


4.1 ± 1.2 


2.0 ±1.5 


0.5 


4.5 


59.0 




0.50- 


-0.65 


41 


3.6 ±1.0 


3.0 ± 1.1 


2.9 ±1.4 


2.3 


2.7 


4.0 




0.65- 


-0.85 


32 


5.3±1.1 


2.1 ± 1.3 


3.2 ±1.6 


0.05 


19.0 


9.0 



Table 2. Mean bin flux densities and K-S probabilities that the Herschel fluxes of objects in luminosity bins are drawn from the background distribution, as a 
function of wavelength. Notes as for Table[TJ 



Catalogued 


Range in 


Objects 


Mean bin flux density (mly) 


K-S probability (%) 


sources? 


logical. 4) 


in bin 


250 


350 fim 


500 fim 


250 fim 


350 fim 


500 fim 


Included 


21.0-23.7 


27 


290.6 ± 1.2 


118.0 ± 1.4 


42.2 ±1.7 


< 10" 3 


< 10" 3 


< 10" 3 




23.7-24.3 


32 


2.4 ± 1.1 


0.5 ± 1.3 


-1.3 ±1.6 


28.3 


69.8 


12.5 




24.3 - 24.6 


35 


11.5 ±1.1 


10.7 ± 1.2 


5.6 ±1.5 


0.010 


0.05 


1.1 




24.6 - 25.0 


40 


9.2 ±1.0 


4.5 ± 1.1 


2.3 ±1.4 


0.8 


10.3 


48.4 




25.0-25.6 


38 


7.1 ± 1.1 


4.7 ±1.2 


4.0 ±1.5 


0.9 


12.2 


5.1 




25.6-27.3 


15 


9.1 ± 1.7 


7.2 ± 1.9 


6.8 ±2.3 


0.10 


0.4 


1.2 


Excluded 


21.0-23.7 


7 


79.2 ±2.5 


42.7 ±2.8 


15.1 ±3.4 


15.2 


39.2 


99.5 




23.7-24.3 


32 


2.4 ±1.1 


0.5 ± 1.3 


-1.3 ±1.6 


28.3 


69.8 


12.5 




24.3 - 24.6 


31 


6.6 ±1.2 


7.6 ± 1.3 


4.7 ±1.6 


0.07 


0.3 


10.8 




24.6 - 25.0 


37 


2.5 ±1.1 


1.8 ± 1.2 


2.0 ±1.5 


5.0 


16.1 


42.9 




25.0-25.6 


35 


3.2±1.1 


2.4 ±1.2 


3.0 ±1.5 


3.2 


32.6 


14.1 




25.6-27.3 


14 


7.7 ±1.7 


5.9 ± 1.9 


6.6 ±2.4 


0.3 


0.8 


1.8 



sities using a maximum-likelihood fit. Serjeant & Hatziminaoglou 
fitted a Gaussian to the flux density distribution of each bin, but 
as Gaussian fits were rather poor to some of our bins, we elected 
to use a lognormal distribution in flux density (which also has the 
merit of being constrained to positive flux density values), applying 
a prior that is uniform in log space to the mean. Rather than assum- 
ing an underlying Gaussian distribution of the noise, we used the 
actual noise distribution in the data (as determined from our ran- 
dom flux measurements) and used Monte-Carlo simulation to de- 
rive the probability distribution of flux densities (which, after the 
addition of noise, can be negative) and thus the model likelihood 
for any given set of parameters. This allowed us to determine the 
best-fitting values for the intrinsic distribution in flux densities us- 
ing a Markov-Chain Monte Carlo (MCMC) algorithm (see Mullin 
& Hardcastle 2009 for a description of the code). Having done this, 
we took noise-weighted means of the luminosities in each bin in the 
manner described by Serjeant & Hatziminaoglou (2009), for sim- 
plicity at this stage making the approximation that both the intrinsic 
scatter in the distribution and the noise were Gaussian so that they 
could be added in quadrature to determine the weights. (We note 
that the overall result of this procedure is reassuringly similar to 
what we obtain if we simply take the mean of the luminosities us- 
ing only the measured noise values; the results are not dominated 



by the weights derived from the intrinsic scatter estimates.) As in 
Section |4~T1 the luminosities were calculated from the measured 
250-fim fluxes (since the detections of stacks are most significant 
in this band, as discussed in Section l4~2t on the assumption of a 
grey-body template SEE0- 

Results are plotted in Fig. [3] Looking first at the results of bin- 
ning by radio luminosity (top-hand panel) we see that, as expected, 
the lowest-luminosity radio bin has a reasonably high mean FIR lu- 
minosity around 10 11 Lq. The highest FIR luminosities of objects 
in this bin are comparable to those of starbursts like Arp 220, which 



2 The ratio of the stacked flux densities for these objects given in TablefT] 
allows us to make a rough check of the temperature of the grey-body model 
assumed in determining the luminosities. For temperatures in the range 20- 
30 K and an assumed = 1.5, the ratio of flux densities at 250 and 350 
fim, where we have the best statistics, should be in the range 1.9-2.7, with 
only a slight dependence on redshift. We see that this is broadly consistent 
with the flux densities in the stacks, although some have ratios closer to 
unity, which would imply even lower temperatures. However, given that we 
have no temperatures for individual sources and so must adopt a single- 
temperature model, it is reassuring that most of the flux density ratios are 
roughly consistent, within their uncertainties, with the T = 26 K model 
that we are using. 



Herschel-ArLASV Radio-selected galaxies 7 



Table 3. Mean bin flux densities of the comparison galaxy population and K-S probabilities that the Herschel fluxes of objects in redshift bins are drawn from 
this galaxy population, as a function of wavelength. 



z range 


Radio-loud 


Mean galaxy flux density (mjy) 


K-S probability (%) 


Comparison 




objects in bin 


250 fim 


350 /im 


500 /urn 


250 fim 


350 (im 


500 ^tm 


galaxies 


0.00-0.10 


15 


22.3 ±0.2 


11.0 ±0.3 


4.7 ±0.3 


< 10~ 3 


< 10 -3 


< 10 -3 


770 


0.10-0.20 


14 


15.5 ± 0.1 


7.4 ±0.1 


3.0 ±0.2 


1.4 


2.0 


4.6 


2913 


0.20 - 0.35 


39 


12.6 ± 0.1 


6.4 ±0.1 


2.9 ± 0.1 


7.3 


4.7 


6.5 


5481 


0.35 - 0.50 


40 


7.9 ±0.1 


4.5 ±0.1 


2.0 ±0.1 


92.5 


65.6 


98.7 


14660 


0.50-0.65 


43 


6.9 ±0.1 


4.4 ±0.1 


2.2 ± 0.1 


41.9 


97.5 


35.1 


11905 


0.65 - 0.85 


36 


9.5 ±0.1 


6.5 ±0.1 


3.1 ± 0.2 


90.6 


58.7 


80.9 


2971 



suggests a comparably high star-formation rate. The second bin has 
a much lower luminosity (which, given that this bin is not detected 
in our K-S test analysis, should be considered to be an upper limit) 
corresponding to star formation rates < 5Mq year -1 . This is the 
effect of moving from a selection criterion that selects mostly star- 
bursts to one that selects mostly AGN, and it immediately shows 
that the lowest-luminosity radio-loud AGN tend to have little or 
no star formation. The remaining radio bins have FIR luminosities 
somewhat higher than the mean in the first bin, and show at most a 
weak positive trend with radio luminosity; these luminosities would 
correspond (using the starburst relation from Kennicutt (1998), and 
bearing in mind the many caveats associated with doing so) to total 
star formation rates between 50 and 100 Mq year -1 , which might 
be associated either with the host galaxy of the radio source itself 
(in which case these would be high star formation rates for quies- 
cent ellipticals) or with nearby galaxies in a host group or cluster. It 
is important to bear in mind that the star-formation rates we quote 
here, and the values plotted in Fig. [3] will be affected by systematic 
uncertainties in the correction from 250-/im rest-frame flux to to- 
tal mid-IR luminosity, as discussed at the start of this Section. The 
relative star-formation rates should be robust, but the absolute nor- 
malizations might be systematically wrong by a significant factor. 
We note, however, that the star-formation rates we derive are quite 
consistent with those derived by Seymour et al. (2010) in a study of 
radio-loud objects in the Herschel Multi-Tiered Extragalactic Sur- 
vey (HerMES). 

The results of luminosity stacking must be interpreted in the 
light of the stacking in redshift bins (Fig. [3] lower panel). Here we 
perform an identical stacking analysis using the normal galaxy pop- 
ulation, selected as defined in the previous Section, as a control. As 
we expected given the results of the K-S tests, we see that the FIR 
luminosities of radio-selected objects are much higher than those 
of comparably selected galaxies in the lowest two redshift bins, 
corresponding to the objects likely to be starburst galaxies. How- 
ever, above this redshift, we see very close agreement between the 
radio-loud objects and the general galaxy population. The differ- 
ence between the two even at the highest redshifts (and therefore 
highest radio luminosities) is little more than la with the current 
data. 

Finally, we can also compare with the properties of optically 
selected quasars from the Bonfield et al. (2010) sample in the same 
redshift bins. There are only 67 unique quasars in this sample in 
the z < 0.85 redshift range, so the sample size is considerably 
smaller than for our radio-selected objects; there are also no objects 
with z < 0.10. However, the individual redshift bins are all signif- 
icantly detected at 99 per cent confidence or better on K-S tests, 
so we can validly use these small samples to compare with the ra- 



dio galaxies. None of the quasars in this redshift range is identified 
with a radio source in our sample. We see (Fig. [3] lower panel) that 
there is a clear trend for the luminosities of these quasars to lie sig- 
nificantly above the luminosities of both the normal galaxies and 
the radio-selected objects, though this tendency appears weaker at 
higher redshift. 



5 DISCUSSION AND CONCLUSIONS 

The most obvious conclusion to be drawn from Section |4~3l is that 
the FIR properties of radio galaxies and those of a comparably- 
selected population of radio-quiet objects are very similar, given 
the data that are presently available to us. The selection of our 
comparison sample could clearly be improved with more data; for 
example, our radio non-detections at high redshift will clearly in- 
clude some lower-luminosity radio-loud objects, while we have not 
attempted to use the available optical data to select exclusively 
elliptical galaxies. But taking the observations at face value, we 
see no evidence that radio galaxy hosts behave any differently in 
the FIR from a matched population of radio-quiet galaxies. We do 
not believe that the incompleteness of our optical identifications of 
the NVSS/FIRST sources should affect this conclusion: the miss- 
ing objects are likely either quasars (and therefore undesirable be- 
cause of the possibility of contamination by non-thermal emission) 
or higher-redshift objects (see below). 

Should we be surprised by our result? We begin by noting that 
we do not expect any effect from any radiatively efficient AGN in 
these radio galaxies — the typical mid-IR luminosities (at around 
15 /im) of radio-loud AGN of comparable radio luminosity to the 
most powerful objects in our sample are 1-2 orders of magnitude 
lower than what we observe in the present sample, if they are de- 
tected at all (Hardcastle et al. 2009) and the emission from the torus 
of radio galaxies believed to be hosting an obscured AGN would 
be expected, and is observed, to peak in the mid-IR (e.g. Haas et 
al 2004). We also do not expect to see any synchrotron contami- 
nation in the FIR, bearing in mind that the typical flux density of 
our sources at 1.4 GHz is a few mjy while the flux density at 250 
/im (1.2 THz) of detected sources is an order of magnitude higher, 
and that, as discussed above, our selection criteria should strongly 
favour steep-spectrum objects. We are thus safe to interpret the FIR 
luminosities as telling us about star formation, subject to the usual 
caveats about young stars being the dominant source of dust heat- 
ing in these objects (Kennicutt 1998). 

The picture appears then to be that the host galaxies of these 
low-luminosity radio-loud AGN have, on average, no more - or 
less - star formation than the general population selected to have 
similar K-band magnitudes. This is the first time it has been possi- 
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Figure 3. Far-IR luminosity, derived from the measured 250- fim luminosity, as a function of (top) 1 .4-GHz luminosity and (bottom) redshift for the radio-loud 
sample, showing both individual objects and the results of stacking as described in the text. Green crosses show identifications in the catalogue of S10; red 
stars show individual objects down to 2<r, and black dots show 2cr upper limits. Error bars on individual luminosities are not plotted for clarity. Blue solid lines 
in both panels are the results of stacking and show the mean luminosity and associated error bars in each of the radio luminosity or redshift bins defined in the 
text. Cyan solid lines in the lower panel show the results of a similar stacking analysis applied to the comparison sample of normal galaxies described in the 
text, and orange solid lines show the results of stacking the quasar sample also described in the text. The right-hand axis shows the conversion between far-IR 
luminosity and star-formation rate, assuming a Kennicutt (1998) starburst relationship; the locus of the limits in the bottom figure thus shows the limiting 
star-formation rate to which we are sensitive for individual objects. 
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ble to investigate this with a large sample in the FIR at these low 
radio luminosities; earlier work studying high-z radio galaxies and 
radio-loud quasars with SCUBA (e.g. Archibald et al.2001; Willott 
et al. 2002; Rawlings et al. 2004) probed radio luminosities which 
were almost all significantly larger than those studied here (for ex- 
ample, all the sources in Archibald et al. have L1.4 > 4 x 10 28 W 
Hz - ), while studies of low-z radio galaxies have generally used 
shorter wavelengths, such as the 60 /im of IRAS (e.g. Yates & Lon- 
gair 1989; Hes et al. 1995) or the 70 /im of Spitzer (e.g. Dicken 
et al. 2009) where the general consensus seems to be that AGN- 
related thermal and non-thermal emission dominates the measured 
flux densities. Our results will therefore be useful for updating the 
latest models of far-infrared emission emanating from AGN and 
their hosts (e.g. Wilman et al. 2010). 

The result above is certainly consistent with our general pic- 
ture of the properties of radio-loud AGN and their hosts. We know 
that at z — the vast majority of low-luminosity radio galaxies are 
hosted by quiescent ellipticals (e.g., Lilly & Longair 1984; Best et 
al. 2005); in these low-power objects there is in general little evi- 
dence for major, gas-rich mergers and so we would expect to see 
little or no ongoing star formation (cf. Kaufmann, Heckman & Best 
2008). But it is interesting to ask how this can be reconciled with 
the fact that there does seem to be a strong association between 
AGN activity and star formation when the AGN are selected at 
other wavebands (e.g. Serjeant et al. 2010). The answer is probably 
related to the two populations of radio-loud AGN, low-excitation 
and high-excitation radio galaxies, described in Section[T] The ra- 
dio luminosity functions of low-excitation and high-excitation ra- 
dio galaxies clearly have different slopes, since low-excitation ob- 
jects are strongly numerically dominant at low radio luminosities 
but almost completely absent at the highest luminosities: the tran- 
sition in terms of numerical dominance seems to take place at ra- 
dio luminosities comparable to, but probably slightly above, the 
conventional Fanaroff-Riley break. Our radio-loud sample is thus 
numerically completely dominated by objects that are likely to be 
low-excitation radio galaxies, which we would expect, based on the 
models discussed in SectionQ] to have little or no association with 
star formation. Only in the highest radio-luminosity bin might we 
expect to see substantial numbers of high-excitation radio galaxies 
with evidence for star formation, and the statistics there are poor. It 
should be noted that Fig. [3] clearly shows that there are individual 
objects with both high radio luminosities and high inferred star for- 
mation rates at all radio luminosities above a few x 10 24 W Hz -1 ; 
these may well be high-excitation radio galaxies associated with 
gas-rich mergers, and they certainly account for the relatively high 
mean FIR luminosities and star-formation rates in the high-radio- 
luminosity bins in that Figure. The idea that high-excitation AGN, 
if present, would be associated with higher star-formation rates is 
borne out, at least qualitatively, by the position of the quasars from 
the sample of Bonfield et al. (2010) in the lower panel of Fig. [3] 
though we caution that we have necessarily made no attempt to 
match these in host galaxy properties to the radio-selected or nor- 
mal galaxies. 

With deeper optical/IR data, better spectroscopy and radio 
data from our GMRT observations, plus broader coverage with 
Herschel (all of which will be available in the near future), we 
can test this qualitative picture in a number of ways. With more 
objects of higher radio luminosity, we will be able to break the 
luminosity-redshift degeneracy inherent to the present sample (the 
results presented here give no new information on whether the far- 
IR luminosity of radio galaxy hosts depends principally on redshift 
or on radio luminosity, a problem first identified by Rawlings et 



al. 2004). The SDP data represent only l/40th of the 570 deg 2 of 
the full H-ATLAS data, so our eventual sample size will increase 
by a large factor: in particular, we expect ~ 6 FRIIs per square 
degree (Wilman et al. 2008) up to the highest redshift so the full 
H-ATLAS will give us > 3000 powerful sources. In addition, even 
with the current data we have optical/IR identifications for less than 
half of the original radio sample (Section [3} so in principle many 
high-redshift radio sources are already there and just await iden- 
tification; the availability of VIKING data (1.4 mag deeper than 
UKIDSS-LAS) will allow us to identify powerful radio sources out 
to z ~ 2.5 (Jarvis et al. 2001; Willott et al. 2003) and so we will 
be able both to extend the stacking study and, probably, to increase 
the number of Herschel-detecled sources, allowing a study of the 
detected population as a function of radio luminosity and morphol- 
ogy. Secondly, the model outlined above implies that we would ex- 
pect any excess FIR luminosity in radio galaxies to be associated 
not only with indicators of star formation at other wavelengths (e.g. 
optical galaxy colours) but also with indicators of high-excitation 
AGN activity such as strong narrow emission lines (which can be 
investigated with GAMA and other spectroscopy), a mid-IR ex- 
cess (which might be visible in survey data from the Wide-field 
Infrared Survey Explorer, WISE) and obscured nuclear X-ray emis- 
sion (many H-ATLAS fields have available X-ray data). With the 
larger sample and the wider availability of multi-wavelength diag- 
nostics provided by the full H-ATLAS project, it should be possible 
to carry out a definitive test of this model. 
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